Electronic exchange in quantum rings 
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Quantum rings can be characterized by a specific radius and ring width. For this rich class 
of physical systems, an accurate approximation for the exchange-hole potential and thus for the 
exchange energy is derived from first principles. Excellent agreement with the exact-exchange results 
is obtained regardless of the ring parameters, total spin, current, or the external magnetic field. 
The description can be applied as a density functional outperforming the commonly used local-spin- 
density approximation, which is here explicitly shown to break down in the quasi-one-dimensional 
limit. The dimensional crossover, which is of extraordinary importance in low-dimensional systems, 
is fully captured by our functional. 
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Ring-shaped quantum systems such as semiconductor 
quantum rings (QRs) have attracted wide interest both 
theoretically and experimentally. Presently, QRs can be 
fabricated using a variety of techniques.— 1^1^ The tunabil- 
ity of size, shape, and electron number in QRs suggests 
applications in the field of quantum-information tech- 
nology. In particular, one can exploit the well-known 
Aharonov-Bohm effect* by external magnetic fields,^ or 
control the electronic states by short laser pulses within 
the decoherence time.^ 

Many-electron QRs have been studied theoretically us- 
ing various approaches, e.g., model Hamiltonians,— ex- 
act diagonalizatioufSiiS quantum Monte CarlO )^°i^^ and 
density- functional theor y -^^'^'^1^* (DFT). Within spin- 
DFT (SDFT) or current-SDFT applied to QRs, the ex- 
change and correlation energies and potentials are com- 
monly calculated using the two-dimensional (2D) lo- 
cal spin-density approximation^^ (LSDA). For relatively 
wide QRs, LSDA performs reasonably well^^ii^ which is 
not surprising in view of the good performance of the 
LSDA in the case of 2D quantum dots.^^ However, in the 
quasi-one-dimcnsional limit the 2D-LSDA is expected to 
fail similarly to the breakdown of the three-dimensional 
LSDA in the quasi-2D limit. Hence, in order to benefit 
from the efficiency of DFT methods in various QRs, accu- 
rate density functionals for exchange and correlation are 
needed. Methods based on exact-exchange (EXX) func- 
tionals seem an attractive possibility.^^ However, most 
computational schemes exploting EXX for finite systems 
suffer from numerical problems p^ii^ which, ultimately, 
prevents the method from being applied to large electron 
numbers. 

In this paper we derive an accurate method to calcu- 
late the exchange-hole potential and the exchange en- 
ergy of QRs. The resulting density functional is sim- 
pler than EXX, but yet compatible with correlation func- 
tionals based on the modeling of the correlation holc.^" 
Our derivation follows the strategy originally proposed 
for atoms by Becke and Roussel,^"'^ in which the averaged 
exchange hole of a suitable single-electron wavefunction 
is adapted to a general TV-electron system by examining 



the short-range behavior of the exchange hole. Recently, 
a similar approach has been used to develop exchange 
functionals for finite 2D systems,^^ which can be repro- 
duced as a special case of this work. In the numerical 
examples we demonstrate the accuracy of the functional 
against exact-exchange results and underline the consid- 
erable improvement over the LSDA, particularly in the 
quasi- ID limit. 

In the framework of SDFT the exact exchange energy 
functional of the spin densities p| (r) and pi (r) can be 
written in (effective) atomic units^* (a.u.) as 



dr2^^^/ij(ri,r2), 



ri -r2 



where 
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is the exchange-hole function. Here we assume that the 
noninteracting ground state is nondegenerate and hence 
takes the form of a single Slater determinant constructed 
from the KS orbitals, V'fc.cr- 

Next we look for an approximation for the cylindrical 
average of the exchange hole. It is defined around ri with 
respect to s = r2 - 

/i^(ri,s) 
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In the following, we compute ft-J exactly for a single- 
electron wavefunction of a QR. We consider a 2D QR 
defined by a radial confining potential of the forni^^ 



V{r) 
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where M > and a > are constants. Note that M ~ Q 
corresponds to a harmonic quantum dot,^^ Regarding QR 
fabricatioUjii'^'^ the confinement given above is realistic: 
On one hand, the electrons cannot enter the center area 
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FIG. 1: (color online) Relation between the parameter set 
{a, M) in the external confining potential [Eq. @] and the 
single-electron ring radius ro and width Ar (see the inset). 
The curves correspond to fixed values of M. 



described by the strongly peaked "antidot" [first term in 
y(r)], and, on the other hand, the edge of the QR is de- 
scribed by a soft, parabohc confinement [second term in 
y(r)]. In this tunable model, both the radius and width 
of the QR can be changed independently by varying M 
and a (see below). 

The eigenfunctions and eigenvalues for a single electron 
confined by V{r) in Eq. ([4]) can be solved analytically.^- 
Setting the radial and angular quantum numbers to zero, 
n = m = 0, yields a wavefunction 



" M - 
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(5) 



We emphasize that this expression, normalized for each 
a, is the general single-electron ground-state wavefunc- 
tion of a QR having a potential minimum at r = ro = 
\fM ja, which corresponds to the ring radius. The width 
of the single-electron QR having the energy Eq — can 
be approximated by Ar = \/2/a. In Fig. [T] we show 
explicitly the relation between the single- electron radius 
and width of the QR (rg , Ar) and the parameters (M, a) 
in the external confining potential given in Eq. 

The exact exchange-hole function for the single- 
electron wavefunction in Eq. ([5]) becomes 



2(A/+1) 
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we first set ri = r and r2 = r + s. As a result we get 
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where l^iz) is the zeroth-order modified Bessel function 
of the first kind, and M is now an integer. 

Next, we apply h% -Y as a genera? model for the averaged 
exchange hole of an iV-electron system.— The model is 
required to locally reproduce the short-range behavior of 
the exchange hole. Therefore, we parametrize o? and r^ 
into functions of r by setting c? a{r) and r^ b{r). 
Equation ([7]) can be now rewritten as 
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where M* — 2(Ar*/ro)^^ is the nearest integer describ- 
ing the characteristic ratio between the effective width 
and the radius of the TV-electron QR. As we will show be- 
low, rg and Ar* can be extracted from the spin density. 
As a consequence, M* can be considered as a (spin) den- 
sity functional M*\p^. The short-range behavior with 
respect to s can be obtained from the first two non- 
vanishing terms in the Taylor expansion in Eq. ([3]): 



K{v,s)=p,{v)^s'C%(v) 



where 



2 Pa 



P<y 



(9) 



(10) 



of the exchange hole in 
given reference point r (argu- 



is the local curvature 
2D^^ around the 
ment omitted). Here the (double of the) kinetic 



energy density, 



and the spin- 



dependent paramagnetic current density, jp,, 



2i S 



k=l 



, depend explic- 



rk,a ~ (V^^,^j ^k,a 

itly on the KS orbitals, and hence implicitly on the spin 
densities. Defining y(r) := a(r)6(r), the zeroth-order 
term in Eq. ^ yields 



Pa 



M'-v 



ttM*\' 

and the second-order term gives 



(11) 



To calculate the cylindrical average as defined in Eq. ^ , 



ttM*\ 
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Combining Eqs. ((TT|) and (fT2l) leads to 



M*\y-^^' +^)[{y-M*f-y\ey = ^, (13) 

from which y can be solved numerically. Now, we can 
compute a and h and thus the averaged exchange hole 
from Eq. ([5]). The exchange- hole potential is given by 



U^{r) = -27r / dshl{a,b\s), 



(14) 



from which the exchange energy can be calculated as 



(15) 



From Eqs. ([5]) and (fTi|) it can be shown that our calcula- 
tion scheme preserves also the exact long-range behavior: 
;7^(r^cx))^-l/r.24 

To summarize our method, the calculation of the ex- 
change energy of an A^-electron QR consists of the fol- 
lowing steps: 

(1) Use the KS orbitals ipk,^ to calculate the spin den- 
sities Pa, the local curvature of the exchange hole 
from Eq. and to estimate M*. 



(2) Compute y{r) numerically from Eq. 

(3) Calculate a{r) from Eq. (fTT|) and b{r) from the re- 
lation b{r) = y{r)/a{r). 

(4) Calculate the averaged exchange-hole function from 
Eq. dHl). 

(5) Finally calculate the exchange-hole potential and 
the exchange energy from Eqs. ([14]) and p^ . re- 
spectively. 

Next we test the performance of our functional in a few 
examples. As the reference results we use the exchange- 
hole potentials and exchange energies of the EXX calcula- 
tions within the rather accurate Krieger-Li-Iafrate (KLI) 
approximation.^^ We compare the results also to the 2D- 
LSDA. Both the EXX and LSDA results are obtained 
using the octopus code;^ The converged EXX orbitals 
are used as the input KS orbitals in our functional (and 
in the LSDA). Alternatively, also the LSDA orbitals can 
be used as input, which in most cases leads to only minor 
changes in the results. 

Figure [5] shows the exchange-hole potentials for two- 
electron singlet states of two QRs defined by (a) 
(M, a)=(l, 0.5) and (b) (M, a) = (9, 3), respectively. Note 
that here we have set M* = M as the first approxima- 
tion. In both cases we find excellent agreement between 
the EXX (dashed lines) and the present work (solid lines). 
Both the large-r limit, where decays as — l/r, and also 
the r limit, are correctly reproduced. 

The exchange energies per particle of the LSDA (dot- 
ted lines in Fig. [J), which are directly comparable to 
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FIG. 2: (color online) Exchange-hole potentials of two- 
electron quantum rings (total spin 5 = 0) defined by param- 
eters (a) M = 1, a = 0.5, and (b) M = 9, a = 3. The dashed 
line shows the exact-exchange (EXX) result, the solid line is 
the result of the present work, and the dotted line corresponds 
to the exchange-energy density in the local spin-density ap- 
proximation (LSDA). 



U^, deviate significantly from the EXX and from our 
approximation. The inability of the LSDA to yield the 
correct shape of the curve is due to the simple density- 
dependent expression of the exchange energy per particle 



in the LSDA, i.e., e 



LSDA 



1/2 

cx pj 



leading to differences 



also in the exchange energy. In the case of Fig. [2Ia), 
we find £;LSDA ^ _o.389, whereas E^^^ = -0.409 and 
^.present ^ _o.408. When the ratio Ar/ro is decreased 
to one third, i.e., M — 1 ^ 9, the deviation of the LSDA 
becomes more pronounced as shown in Fig. [^b). Now 
we find £:LSDA ^ _i.502 vs. E^^^ = £;p>-°^™t = -1.300. 
In other words, the above change in M increases the rela- 
tive error of the LSDA exchange energy from 5% to 16%. 
Hence, these results demonstrate the breakdown of the 
2D-LSDA in the quasi-one-dimensional limit. 

At this point we make two remarks of the numerical 
procedure. First, Eq. (|13p has two solutions for y(r), 
from which we choose the smaller one for r < r[p,nax] 
and the larger one otherwise. Since these two solutions 
do not generally coincide at r = r-[pniax], we extrap- 
olate U^{r) around this point as visualized in the in- 
set of Fig. [51 This smoothing procedure is done such 
that Ex is not altered. Second, the application of the 
overall scheme might be numerically cumbersome in the 
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FIG. 3: (color online) Same as in Fig. [2] but for (a) a fully 
spin-polarized current-carrying six-electron quantum ring at 
B = 6 T and for (b) a zero-current 12-eIectron quantum ring. 



polarized current-carrying six-electron QR at B = 6 T 
and for (b) a 12-electron QR. Similarly to the previous 
examples, we find excellent agreement with the EXX. 
The six-electron case yields E^^^ = -2.226, EP"''""^^ = 
-2.238, and E^^^^ = -2.114. In the case of 12 elec- 
trons, we plot results for both M* = Af = 3 and for 
M*[p] = 2 « 2(Ar*/rQ)~^, where the effective radius 
corresponds to the point where the cumulative density 
reaches 50%, and the effective width Ar*, centered at 
Tq , covers 90 % of the total density. We find that in the 
latter case the agreement with the EXX is better. 

To conclude, we have derived, from first principles, an 
accurate and general approximation for the exchange- 
hole potential and hence for the exchange energy in quan- 
tum rings. Excellent agreement with the exact-exchange 
results is obtained regardless of the ring geometry, num- 
ber of electrons, spin polarization, and currents. More- 
over, we have demonstrated that, in contrast to the local- 
density approximation, our functional can deal with the 
physically relevant dimensional crossover between two 
and one dimensions. Our approach is suitable for the 
development of correlation functionals by considering the 
exact properties of the corresponding correlation-hole po- 
tentials3 
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small-r regime when M (and/or M*) is large and the 
density is very small. To overcome this problem, we 
take the advantage of the exact property of our model, 
U^{r ^ 0) = -(2M)! (M!)-2 2-2*^ + 0{r^), which 
is valid for a general A^-electron case with M > 1^ 

To demonstrate the generality of our method, we show 
in Fig. [3] the exchange-hole potentials for (a) a spin- 
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